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ABSTRACT 

We present an analytical model of a magnetar as a high density magnetized 
quark bag. The effect of strong magnetic fields (B > 5 x 10 16 G) in the equation 
of state is considered. An analytic expression for the mass-radius relationship is 
found from the energy variational principle in general relativity. Our results are 
compared with observational evidences of possible quark and/or hybrid stars. 

Subject headings: dense matter-magnetic fields-methods: analytical-stars: indi- 
vidual(Quark Stars)-stars: magnetars 



1. INTRODUCTION 

The fundamental aspects of the physics involved in the des cription of the matter inside 



a white dwarf are well understood ( IShapiro &: Teukolskylll983l ). but in the case of neutron 



stars the situation is rather different because the equation of state (EoS) of neutron matter 
at very high densities is still unknown. 

The interior of a neutron star is an astrophysical laboratory in which matter is com- 
pressed to high densities. The compression of matter several times the saturation nuclear 
matter density, po, may produce a phase transition from nuclear to quark matter, i.e., an 
unconfined quark - gluon plasma. In addition, under suitable circumstances, a conversion 
d — > s quarks may happen through weak interactions, leading to what has been called 
strange quark matter (SQ M). It has been stated that SQM may be the absolute ground 



state of strong interactions ( IBodmerlll97ll : IWittenlll984j ). although such hypothesis has not 



been confirmed yet. The natural scenario where SQM could occur is the inner core of neutron 
stars. Hence, if the SQM hypothesis is true, some neutron stars could be either hybrid stars, 
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which have quark cores surrounded by a hadronic shell, or quark stars. Al ready 40 ye ars 
ago, the existence of quark stars in hydrostatic equilibrium was suggested by lltohl (ll970l ) in 
a preliminary work. 

On the other hand, it is well known that at the surface of neutron stars there exist 
magnetic fields of the order of 10 12 — 10 13 G. Compact stars with ultra strong magnetic fields 
(10 2 — 10 3 larger than those of a typical neutron star) are call ed magnetars. In s uch objects 
the magnetic field at the surface could be higher than 10 15 G (jde la Incerall2009l ). 



The knowledge of the magnetars composition would help explain some astrophysical 
phenomena. Soft gamma - ray repeaters (SGRs) and anomalous X - ray pulsars (AXPs) 
have been interpreted a s evidences of magnetars. Besides some authors (jCheng fc Daib 



2002 



Ouyed et al.l 120071 ) claim that magnetized hybrid or quark stars could be the real 
sources of SGRs and AXPs. The M-R relationship tells us how matter composing the 
star behaves under compression, providing information about its composition. Several EoS 
for neutron, hybrid and quark stars have been proposed but none of them is co nclusive 
rtDouchin fc Haensell boOll ; lLattimer fc Prakash I boOll . 12007k Ipzel fc Psaltisl hood ). Each 
EoS produces a different mass-radius (M-R) relationship which can be contrasted with the 
available observational data in order to test its range of validity and/or set bounds on some 
parameters. At this particular point astrophysical studies become of great importance since 
they could shed some light in understanding fundamental aspects of matter: microphysics 
could be inferred from macrophysics. Here lies the great importance of the studies related to 
ultra-compact objects. For instance lLattimer &: Prakash I (120011 . 120071 ) contrast some M-R 
relationships obtained theoretically for different EoS. Varying some parameters a difference 
of 4% — 10% and 10% — 15% in determining the maximum radius R ma x and mass M max , 
respectively, is shown for the same EoS. 



Several papers (jChakrabarty fc Sahulll996l ; lGonzalez Felipe fc Perez Martmezll2009l : lPerez Martinez et 



2010) study the M-R relationship of highly magnetized quark stars (HMQS) through numer- 
ical integration of the Tolman-Oppenheimer-Volkoff (TOV) equation for different EoS. Al- 
though most studies of quark stars properties have used such method, iBanerjee et. all (120001 ) 
have obtained a maximum mass and radius for unmagnetized quark stars analytically by 
using a non-relativistic gravitational treatment. 

Approximate analytical solutions play an important role in the astrophysical analysis, 
giving a keener insight than the numerical solutions. Moreover, they may be used as a 
testing point to check if the numerical scheme is accurate and also they are the first step in 
the comparison between theory and observation. Indeed, an approximate analytical solution 
for M-R relationship may be all that is required when comparisons with observational limits 
that determine the confidence contour for the mass and radius are performed. Besides, in 
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the high - density EoS the uncertainties are of the same order or larger than the errors in 
the variational method. 

The appropriate treatment for quark stars should be relativistic, since the existence of 
a maximum mass is assoc iated with the b ehavior of a relativistic gas and general relativistic 
corrections are dominant (IWeinberglll972l ) . In this paper, we shall use the general relativistic 
energy variational principle described by iNaurenberg fc Chaplind ( 119731 ) to obtain an ana- 
lytic approximate formula for the mass, radius, and baryonic number of an HMQS. Quark 
stars are particularly suitable for a variational treatment since their density profile resembles 
a constant mass density star. We shall model an HMQS assuming quark matter within high 
density regime in the framework of a modified MIT Bag model EoS. We also assume that 
the magnetic field B is low enough to be treated like a correction in the EoS (B << /i 2 , with 
/j, being the baryon chemical potential) although, as we will see in the following sections, this 
is not a strong restriction. 

The paper is organized as follow. In Section [2] we calculate the thermodynamical quan- 
tities of the system and analyze the stability of quark matter with respect to decomposition 
in baryons. In Section [3] we provide the analytic relativistic M-R relationship and compare 
our results with the observational data. We also check the dynamic stability of the star by 
calculating the adiabatic index and the speed of sound . In Section HI we present a summary 
of our main results and conclusions. 



2. High density quark matter within a strong magnetic field 

In this section we shall discuss the analytic approximations to the SQM EoS in the 
presence of an uniform magnetic field B \\ z. Within the framework of the MIT Bag model, 
we assume three massless quarks u, d, and s, neglecting mediated interactions between them. 
We also consider that the strong magnetic field is a small contribution to the total energy, 
a fact that will be checked later. 



2.1. Quark Matter in a Magnetic Field 

Let us compute the grand canonical thermodynamic potential Q in the high density 
regime. Due the Landau quantization the phase space volume integral in the momentum 
space is replaced by 



(2vr 



^3 / d'PfiP) = 7^)3 / dPzd V f(p) = g X> - **) /_ 



dp z f(v,Pz), (1) 
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where (2 — 5 u o) means that the zeroth Landau level is singly degenerate, whereas all other 
states are doubly degenerate. The grand canonical potential for each quark in the presence 
of a strong magnetic field is given by 



v=0 



/i V > - - 2^5 In 

V2z/gii? 



(2) 



where ^ = 2 x 3, taking into account spin and color degeneracy, and qi is the absolute value 
of the charge of the particle, q u = 2|e|/3 and qd = q s = |e|/3, with e being the value of 
electronic charge. 

For simplicity, we consider the quark masses m q = 0, which implies that the electrons are 
not present and quarks chemical potential are, as a consequence of equilibrium conditions, 
all equal, /i u = fid = fi s = fi. 

By imposing that 

pi = fi 2 - 2uq t B > 0, (3) 
we can determine the upper limit of the sum v max from 

V ~ = IW (4) 



The series, Equation (j2j), can be approximated using the Euler-MacLaurin formula 



£/0') 

j=0 



/(x)dx + -[/(n) + /(0)] 



- /'(«)] - ^LT(n) - f"(o)] + ^ 



(5) 



where the rem ainder term, R, usually is expressed in terms of periodic Bernoulli polynomials 
(jSpiveyl 120061 ) . and can be estimated by using 



\R\ < 



204) 

(27T) 4 



f IV (u)\ du, 



(6) 



where the Riemann Zeta function £(4) ~ 1.0823. To avoid divergences appearing in the 
third term of Equation (jSJ), in the limit of high densities or negligible quark masses, we apply 
the Euler-MacLaurin formula in the form 

lax 1 



1 



Vti{v)dv + - [Vti{v max - 1) + fii(l)] 



+ 



12 



dv 



I {yraax 1) 



dv 



l(i) 



+ R, 



(7) 



-5 - 



where 



R 



1 

720 



d 3 {li 



I {y-max 1) 



+ R. 



We have considered Equation (JSJ as the remainder term, which gives | R |< 3%. In the limit 
/i 2 ^> 2qiB, the thermodynamical potential can be calculated performing first the integral 
in Equation ([Tj) and then expanding in power series of B. The result is 



a 



4+^ ! 6og|4-3)H-0(B 5 ' 2 ). 



47T 



S7T 



'2/i 2 



The particle density n.j = — ^ is 



n 4 



7T 2 47T 2 /i 



(9) 



(10) 



Note that when B = in Equations (Q and (ITUj) . we recover the usual expressions for a 
non-interacting massless quark gas at zero temperature and zero magnetic field. 



2.2. Equation of State 

With the above results, one can form the modified EoS of SQM in the MIT Bag model. 
Within this framework, the difference between the energy density of the perturbative and 
non-perturbative QCD vacuum is taken into account by the "bag constant" B. Considering 
h = c = 1 we can find the conversion factor between high energy density units and magnetic 
energy density units. We can write 1 MeV ~ 1.6 x 10~ 6 erg and 1 MeV ~ (2 x 
10 cm y 1 , where 1 MeV 4 - 2x 10 26 erg cm . Relating this quantity with the magnetic 
energy density B 2 /8n , the conversion factor for the magnetic field is given by 1.4 x 10 13 G 
I MeV 2 . 

The charge neutrality condition 

2n u = n d + n s (11) 

and the /9-equilibrium condition 

Hu = V>d = Us = yU ( 12 ) 

are automatically satisfied. 

Combining the results of Section [2] we obtain 

n = £ n, + B eff = + J?l (io g¥ £— - s) + Bcff + o(b^), (is) 

i=u,d,s 
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where B fj = |^ + B . Replacing Equation ( II ip in the baryon number density condition, 
n B = \ J2i=u,d,s n * ' we obtain 



- ■ 0i " : '^- ( 14 ) 



Since we work in the T = limit, the energy density is given by 



p = « + aw = % + ^ (i + log^,) + B cff + o(b^), (is) 



whereas the pressure reads 



Note that we are not considering the anisotropy of pressures (jGonzalez Felipe et al.l 120081 ) 
because we are working in the limit of weak magnetic field, /z 2 3> 2qiB. The relation between 
the total energy density, Equation (|15p . and the total pressure, Equation f fl6|) . determines 
the EoS of the system as 

p = 3P + 4Bcff - * ( 2 - Log^) + 0(flV»). (17) 



2.3. Stability Analysis: Strong Interactions 



It is well known that SQM may be stable with respect to decay into nucleons at zero 
pressure and zero te mperature if its ene rgy per baryon ^ is less than the energy per baryon 
of 56 Fe = 930 MeV (IFarhi fc Jaffdl 19841 ). The presence of a magnetic field changes somewhat 
this stability condition. 



At P 
method as 



we can estimate the chemical potential through successive approximation 

.1/4 



KB, B) 



47r 2 B cfr B 1 

— H 

3 3 2 



log- 



B 



(1* 



2 4 / 3 7r\/3B 

which will be replaced in equations (TMl JT7|) to evaluate Contrary to previou s results 
( lAnand fc Singh! Il999t IChakrabartyl Il999i ; iGonzalez Felipe fc Perez Martinez! |2009|) we find 
that the energy per baryon increases with B (equation ( 1TT1) ). The condition ^ < 930 MeV 
is satisfied for magnetic fields B < 4.4 x 10 18 G . However for the stability of the system, 
not only it is necessary to consider the energy per baryon, but also the influence of magnetic 
energy density. We obtain B 2 /8n ~ B for 85 MeV fin" 3 < B < 90 MeV far 3 and 
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B 2 /8tt > B for 57 MeVfm" 3 < B < 80 MeVfm" 3 . Although it is known that the binding 
of the quark stars is provided not by gravitation, but rather by the strong interactions, 
the inclusion of a magnetic field B adds an additional constraint to the stability condition 
through the magnetic energy density and the magnetic pressure. When the latter is not lower 
than B the magnetic field becomes dynamically important. Furthermore if the magnetic 
pressure is of the same order of magnitude than the matter pressure, spherical deformation 
effects should be considered. In addition sufficiently strong magnetic fields can generate an 
anisotropic pressu re distribution inside the HMQS modifying the EOS and consequently the 
M-R relationship (IPaulucci et al.ll201ll ). Therefore for magnetic fields large enough, B e g- will 
be greater than the kinetic energy of the quarks thereby destabilizing the star. Thus we also 
consider the relationship 

(B 2 /8n) B- 1 < 0.1 (19) 

to guarantee both the perturbative treatment of the system and the stability of the star. 
Table [1] shows that the variation of baryon density is quite small for P = when increasing 
the magnetic field from up to B max . 



3. MASS-RADIUS RELATIONSHIP BY VARIATIONAL METHOD 



The energ y variational meth od in general relativity is explained in detail in ([Harrison et al 



(119651 ): see also lWeinbergJ (119721 )). Starting from an uniform density configuration in a spher- 
ically symmetric distribution the total mass M, the baryon number Nb and the radius R of 
the star are given by 



M 

N B 
R 



4 , 
^ P R 3 , 

2Trn B a 3 (x ~ sin x cos x), 
asinx 



(20) 



where p is the mass-energy density and the angle x comes from substituting r = a sin x, 
where a = [3/(87rp) ] 1//2 is the curvature radius in the metric inside the star which adopts 
the following form for the 3-geometry: 



[dx 



2 + sin 2 x (d8 2 + sin 2 



6d(j) 2 )] . 



(21) 



Note that we are using h = c = G = 1. The configuration of maximum density is achieved 
when x — 7r /2. Observe that sin 2 x = 2M/R, x ~ corresponds to the Newtonian 



limit while x = 1T /2 corresponds t o the Schwarzschild one 



configuration has been justified by lNaurenberg fc Chaplind (119731 ) . while IWeinbergl ( 1 19721 ) has 



he u se of a constant energy t rial 
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applied it to white dwarfs. In this latter case, a fair approximation to the Chandrasekhar 
mass was obtained. 

To obtain the equilibrium condition is appropriate to treat x as an independent variable. 
Imposing d M / d\ — for fixed Ng, the equilibrium condition reads 

w = - = C(x), (22) 

where p and P are given by Equations (|T5j) and (fIE|) . and ((x) is a function independent of 
the EoS 3 

(( x ) = 3 cosx ( - cosx S -^-^ ] -1. (23) 

We get an approximate value of C(x) using a Taylor series, (t , around x = . Truncation 
at eighth order gives 

1 2 113 4 1747 fi 689,687 8 

/• = y 2 J y 4 J y — I ! y 8 . (24) 

s 10 2100 63,000 48,510,000 v ; 

The Pade approximant of order (4 , 4) gives a representation of this function that is also an 
approximate analytic continuation beyond the circle of convergence. Thus, C{x) is given as 
a ratio of two polynomials as 

( 23 4 1 ,\/ 5123 2 3002 A -1 
°> = [-6237 X + V0 X ){ 1 -89T0 X + 93T555 X J " (25) 

Imposing (p = w we obtain the only physical solution for x, always positive and fulfilling 
the condition lim^ _> X = , given by 



^3 J (35861 w + 6237 - V786718681 w 2 + 389949714 w + 38900169) 

X = . (26) 

2 V3002 w + 345 

Hence we get an analytical expression of the HMQS mass and radius as a function of the 
baryonic chemical potential. This allows us to obtain the M-R relationship for different Bag 
constant and magnetic field values. In particular in Figured] we show the B = 57 MeVfm -3 
and B = 90 MeVfm -3 cases with and without magnetic field. Note that for the first one, 
when B = 0, the M-R curve coincides with the hadronic star zone. This result can be 
attributed to the fact that the value B = 57 MeVfm" 3 could be too low f or a quark star 



which has been modeled by using the MIT bag model (IZdunik et al. II2000I ). Furthermore, 
in Table |2] we present the results for other values of B. Note that although B max is a typical 
value for a magnetar it slightly decreases the M max if compared with the zero magnetic field 
case. The same result is obtained for R mnT . 
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3.1. Dynamical Stability 



In our model the condition for stable equilibrium is given by d 2 M / d 2 x > 0. For a 
given EoS, it is possible to determine the quark densities and pressures where quark stars 
are stable against gravitational collapse from the condition 

r>r c , 

where the adiabatic index for SQM, T, is given by: 

■ 4 25V (7-6Log w B 5 



riB dP 
P driB 



'2V33, 



3/i 4 - 4tt 2 B 
88 B 2 B rff 7T 2 
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(3/i 4 -47r 2 B cfr ) 2 
+ 0(B 4 ), 



27 (3/i 4 -47r 2 B cff ) 2 
and the critical adiabatic index, T c , for a cold star in general relativity is 



r c = (i + w) 



i + 



(3w + 1; 



(W + 1) n 

■ tan 2 x ~ 1 



6w 



(27) 



(2f 



To get dynamical stability the condition r > | must be satisfied. The intersection between 
T and T c determines when quark star becomes gravit at tonally unstable. We found r c ~ 2.3 
and w c ~ 0.17 for the maximum mass value of the star (Table [2]). 

Another quantity that is related with the stability of the star is the speed of sound c s . 
To satisfy the causality of quark matter, 



(29) 



At very high densities particles become relativistic and the speed of sound should be lower, 
more precisely of the order of 1/ y 7 ^, the speed of sound of relativistic fluids. We found 



VS 27V3fi 4 



+ 0{B 2 ). 



(30) 



This quantity tell us how stiff is the EoS providing information about the compressibility of 
the fluid. The stiffness of the EoS increases when c s is closer to 1. 



4. SUMMARY AND CONCLUSIONS 



In this article we have furnished an analytical treatment to study an HMQS in the 
framework of the MIT Bag model. We have analyzed the stability of quark matter with 
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respect to strong interactions and found a restriction in the stability condition: there is a 
maximum value for the magnetic field, B max , beyond which quark matter becomes unstable. 
In the limit of "weak" magnetic field that we have studied, quark magnetic moments are 
aligned in the same direction of the field and this situation leads to such restriction. This 
result could mean that if the magnetic field strength exceeds that critical value, then quark 
or hybrid stars should not be considered as magnetars. In a more general case, when quark 
masses are take into account, electrons should be considered. For magnetic fields much 
stronger than B max , Landau levels for electrons will increase the energy per particle. However 
this probably will not contribute to modify the iron-condition , 56 Fe = 930 MeV, because 
electron fraction in SQM is already very low. 

We have also found an analytical approximate solution for the M-R relationship. Even 
tough we used very simple physics, our results are in good agreement with the confidence 
contours of available observational data. 

It is important to note that although the uniform energy density regime is a good 
approximation for quark stars, deviations in the determination of M-R relationship may 
occur because in the limit of high densities such approximation is no longer valid. 

Finally, we calculate the adiabatic index and the speed of sound. The critical value for 
th e adiabatic index, which corres ponds to the collapse of the star, is in agreement with that 



of iNaurenberg &: Chaplind ( 119771 ). a pioneering work about quark stars. On the other hand, 



the speed of sound is consistent with the expected values for quark stars. 
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Fig. 1. — Mass- Radius relationship with and without magnetic field for B = 57MeVfm -3 
(solid and dashed line) and B = 90MeVfm -3 (solid and dash-dotted line). The rectan- 
gle with diagonal pattern corresponds to EXO 0748-676, which has been interpreted as 
a hadronic star. Recta ngles with crossed, vert ical, and horizontal patterns correspond to 
quarks or hybrid stars (D rago fc Lavagndl2010l ). The polygon could be a low-mass strange 
star as suggested in (jZhang et al. l2007h . 
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Table 1: Bag Constant, Baryon Density and Magnetic Field Upper Limit to Preserve Quark 
Matter Stability Condition. 



B[MeV/m 



-31 



57 
60 
75 
80 
85 
90 



n B /n 
1.73 ±0.01 
1.80 ±0.01 
2.14 ±0.01 
2.24 ±0.02 
2.34 ±0.02 
2.45 ±0.02 



■Bmax [C] 

4.4 x 10 17 
4.8 x 10 17 
5.3 x 10 17 

5.5 x 10 17 

5.7 x 10 17 

5.8 x 10 17 
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Table 2: Maximum Mass, Maximum Radius and Baryonic Number for Different Bag Con- 
stants. 



B [MeVfm- 3 


B[G] 


7~i r 7 "1 

Rmax [km\ 


M max /M Q 


N B /N & 


57 







12. fO 


2.55 


4.30 


4.4 


x 10 1 ' 


12.06 


2.45 


4.09 


60 







11. oU 


/1Q 


A 1 A 
4. 14 


4.8 


x fO 17 


11.76 


2.48 


4.14 


75 







10.55 


2.22 


3.50 


5.3 


x fO 17 


10.52 


2.22 


3.50 


80 







10.22 


2.15 


3.34 


5.5 


x fO 17 


10.18 


2.15 


3.34 


85 







9.91 


2.09 


3.19 


5.7 


x fO 17 


9.89 


2.09 


3.19 


90 







9.63 


2.03 


3.05 


5.8 


x fO 17 


9.60 


2.03 


3.05 



